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The formation of spacetime singularities is a quite common phenomenon in General Relativity and 
it is regulated by specific theorems. It is widely believed that spacetime singularities do not exist in 
Nature, but that they represent a limitation of the classical theory. While we do not yet have any 
solid theory of quantum gravity, toy models of black hole solutions without singularities have been 
proposed. So far, there are only non-rotating regular black holes in the literature. These metrics can 
be hardly tested by astrophysical observations, as the black hole spin plays a fundamental role in 
any astrophysical process. In this letter, we apply the Newman-Janis algorithm to the Hayward and 
to the Bardeen black hole metrics. In both cases, we obtain a family of rotating solutions. Every 
solution corresponds to a different matter configuration. Each family has one solution with special 
properties, which can be written in Kerr-like form in Boyer-Lindquist coordinates. These special 
solutions are of Petrov type D, they are singularity free, but they violate the weak energy condition 
for a non-vanishing spin and their curvature invariants have different values at r = depending on 
the way one approaches the origin. We propose a natural prescription to have rotating solutions 
with a minimal violation of the weak energy condition and without the questionable property of the 
curvature invariants at the origin. 



I. INTRODUCTION 



Under the main assumptions of the validity of the strong energy condition and of the existence of global hyperboHcity, 
in General Relativity collapsing matter forms spacetime singularities At a singularity, predictability is lost and 
standard physics breaks down. In analogy with the appearance of divergent quantities in other classical theories, 
it is widely beheved that even spacetime singularities are a symptom of the limitations of General Relativity and 
that they must be solved in a theory of quantum gravity. While quantum gravity effects are traditionally thought to 
show up at the Planck scale, Lpi ~ 10"'^'^ cm, making experimental and observational tests likely impossible, more 
recent studies have put forward a different idea [21 [3] : ipi would be the quantum gravity scale for a system of a few 
particles, while the quantum gravity scale for systems with many constituents would be its gravitational radius. In 
these frameworks, even astrophysical black holes (BHs) of tens or millions Solar masses may be intrinsically quantum 
objects, macroscopically different from the Kerr BHs predicted in General Relativity. 

While we do not yet have any mature and reliable candidate for a quantum theory of gravity, more phenomenological 
approaches have tried to somehow solve these singularities and study possible implications. In this context, an 
important line of research is represented by the work on the so-called regular BH solutions 4}|7]. These spacetimes 
have an event horizon and no pathological features like singularities or regions with closed timelike curves. Of course, 
their metric is not a solution of Einstein's vacuum equations, but they can be introduced either with some exotic 
field, usually some form of non-linear electrodynamics, or modifications to gravity. They can avoid the singularity 
theorems because they meet the weak energy condition, but not the strong one. 

The purpose of the present letter is to construct rotating regular BH solutions. This is a necessary step to test these 
metrics with astrophysical observations [5HTU] . The spin enters as the current-dipole moment of the gravitational field 
of a compact object and it is thus the leading order correction to the mass-monopole term. It is not possible to 
constrain deviations from classical predictions without an independent estimate of the spin. However, exact rotating 
BH solutions different from the classical Kerr-Newman metric are very hard to find. In most cases, including all 
the regular BH metrics currently available in the literature, we know only the non-rotating solution. In a few cases, 
we have an approximated solution valid in the slow- rotation limit which is also not very useful for tests. A 
rotating solution in the Einstein-Gauss-Bonnet-dilaton gravity has been recently found numerically in Rcf. [12 , while 
proposals for some rotating quantum BHs have been suggested in [131 HI! ■ 
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II. NEWMAN-JANIS ALGORITHM 



Our strategy is to use the Newman- Janis transformation [15] (for more details, see Ref. [Hj). Roughly speaking, 
the algorithm starts with a non-rotating spacetime and, at the end of the procedure, the spacetimc has an asymptotic 
notion of angular momentum. The starting point is a spherically symmetric spacetime 

ds^ = f(r)dt^ - — - h(r) (de'^ + sin^ Odcf)'^) . (1) 
f{r) 

The first step of the algorithm is a transformation to get null coordinates {u, r, 6, 0}, where 

du = dt-dr/f{r). (2) 
The second step is to find a null tetrad Zj^ = {l^ , n^, m^, fh^) for the inverse matrix in null coordinates 

where the tetrad vectors satisfy the relations 

= m^m^^ = n^n^ = Z^m^ = n^m^ = , ^^n^ = -m^rh^ = 1 , (4) 
and X is the complex conjugate of the general quantity x. One finds 

The f/iird step of the procedure is the combination of two operations. A complex transformation in the r ~ u plane as 
follows 

r-^r' = r + ia cos 9 , m— > u' ~ u — ia cos , (6) 
together with a complexification of the functions /(r) and h(r) of the metric. The new tetrad vectors are 

l'^ = Si^ , n'^ = 6^ - m'^ = — ^ L smOiS^ - 5'^) + 5^, + ^5^] , (7) 

J2h{r') ^ 

where /(r') and h{r') are real functions on the complex domain. This step of the procedure is in principle completely 
arbitrary. In fact, in the original paper, Newman and Janis could not give a true explanation of the procedure if not 
that it works for the Kerr metric with a particular choice of the complexifications. The situation improved with Drake 
and Szekeres in [16], where the authors proved that the only Petrov D spacetime generated by the Newman- Janis 
algorithm with a vanishing Ricci scalar is the Kerr- Newman solution. Using the new tetrad in Eq. (|3|, we find the 
new inverse metric and then the metric. The non-vanishing coefficients of g^^ are 



gun = f{r,0), g„r = gr„ = l, gu^f, = g^u = asin'^ 9 (l - f{r, 9)j , 



g^r = a sin^ 9 , gee 



h{r, 9) , - - sin^ 9 \h{r, 9) + a" sin^ (2 - /(r, 0)) ] . (8) 



The fourth and last step of the algorithm is a change of coordinates. In some cases, we can write the metric in 
the Boyer-Lindquist form, in which the only non- vanishing off-diagonal term is gtrj,. This requires a coordinate 
transformation of the form 

du = dt' + F{r)dr , del) = dcj)' + G{r)dr , (9) 

where 

F[r)^ , /^M) + a^Bin^g _^ ^ 

f{r,9)h(r,9) + a^sm^9 f{r,9)h{r, 9) + sm^ 9 

This transformation is possible only when F and G depend on the coordinate r only. In general, however, the 
expressions on the right hand sides of ( 10 ) depend also on 0, and we cannot perform a global transformation of the 
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form ([9]). If the transformation (|9| is allowed and we go to Boyer-Lindquist coordinates, the non- vanishing metric 
coefficients of the rotating BH metric are: 



9tt = f{r, 9) , gt^ = g^t = a sin^ 9 (l - /(r, 9) 



Hr,9) 



h{r, 9)f{r, 9) + a'^ sin^ 



gee = -h{r,{ 



— sm 



h{r, 9) + a^ sm^ 9 [2 - f{r,9) 



(11) 



In the case of the Schwarzschild solution, we have f{r) — 1 — 2M/r and h{r) — . In the Newman-Janis algorithm, 
we have to choose a complexification of the 1/r and of the term. In general, this prescription is not unique. 
However, since wc know what the Kerr solution is, we know that if we take the following complexification 



r 2\r' 



then this trick works well. The functions /(r) and h{r) become 



f{r)^ }{t,9) = 1- 



2Mr 



h{r) h{r, 6*) = E, 



(12) 



(13) 



where S = + cos^ 9. In this case, the functions F and G in Eq. ( 10 1 depend on r only and we find the Kerr 
solution in Boyer-Lindquist coordinates 



2Mr\ 



4aAfrsin^0 , ,, S , , , 
dt^ H dtdch - —dr^ - Y.d9^ - sm^ 1 



2 I 2 
r + a 



2a^Mr &in^ 9 



(14) 



where A = - 2Mr + a^. 



III. HAYWARD BLACK HOLE 



As first example of regular black hole, we consider the Hayward metric, whose analytic expression is quite simple [5]. 
The line element is given by Eq. ([I]), with the following /(r) 



/(r) = 1 - — , m = m(r) = , (15) 

7'-^ + g"^ 



where M is the BH mass and g is some real positive constant measuring the deviations from the classical Kerr metric. 
Let us note that m(r) may be interpreted as the mass inside the sphere of radius r and approaches AI as r goes to 
infinity. This spacetime is everywhere regular, as can be verified by its curvature invariants: 



j^^^l^llS^l^^ limi?=-?l^, (16) 



72M2g6 (5^6 _ 333 _^ 296) 144Af2 



i?^,^,i?^-P- = ^ ' J i , lim i^^.p.i?'^-"- = . (18) 



The weak energy condition is also satisfied. 

When we apply the Newman-Janis algorithm to get a rotating solution, the key-point is the complexification of 
/(r), as the one for h{r) must be the same of Eq. (13). With the only requirement to recover the Kerr metric for 
g = 0, all the possible complcxifications have the form 

/(r,0) = l-^, (19) 

where, in general, rh = rha^(j{r, 9) is a function of both r and 9 and the complexification is characterized by the two 
real numbers a and /3: 



j,3-|-Q5]-a/2 



- ^ ^3+.S-./2 + „3,,S-,/2 ■ (20) 
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At this point, we can distinguish two classes of solutions. The first class has only the case a = /? = (complexifi- 
cation of type-I): we complexify the 1/r term as in Schwarzschild, without altering the mass term m(r). With this 
choice 

fj{r,e)^l-^-^, (21) 



F and G depend on the coordinate r only, and the final result is the line element ( 14) with rn(r) in Eq. ( 15 1 replacing 
M. The spacetime is of Petrov type D, as we can verify by the presence of the Carter constant for the motion of a free 
particle. The new solution is also everywhere regular for 7^ 0, as can be seen from the expression of its curvature 
scalar, Ricci square, and Kretschmann invariant. We do not report the analytic forms here, but they can be quickly 
obtained with Mathematica and any good package for tensor calculations. At the origin, they reduce to 

24Af 



lim f lim i?") = lim ( lim i?)=0, lim ( lini R] = —, (22) 



lim ( lim RnuR^") = lim ( lim Rn^R^'A = , hm ( lim Rni.R^'"^ = , (23) 

e^any Vr^O j r^O \e^B^-^l2 ^ j r^O \e^7r/2 J 9^ 

Yim {\imR^,p„R^-pA=lim( lim R^^^p^R^'^'A = , lim ( lim R^,p,R^-P''\ ^'^-^ . (24) 

The fact that these curvature invariants assume two different values for r = 0, depending on the way one approaches 
the origin, is a signature of the "de Sitter belt" , absent in the non-rotating metric. It has been already found in some 
non-commutative geometry inspired BHs [14] . 

The weak energy condition, satisfied in the non-rotating case, is violated for a ^ 0. To check it, we can choose 
an orthonormal basis in which the stress-energy tensor is diagonal, T^^^^^' = diag(p. Pi, P2, ^a)- The weak energy 
condition requires p > and p -I- Pi > (z = 1, 2, 3) [T]. The one forms of the dual basis of the orthonormal tetrad of 
the standard locally non-rotating frame are |17j 



e(0) ^ 



5t0 

9tt 

94>4> 



1/2 

dt, eW = |.9„|i/2dr, e(2) = (-5,,)i/2d0, e^^) = -^^dt + (-.g,,)^/^ d</) . (25) 



In this frame, the Einstein tensor has an off-diagonal element G'*^°^'^^\ so we need a transformation (e(°\ e*^'^-') — > 
(e'(°\ e'*^'^'') to get a diagonal tensor. At this point, the density p is the eigenvalue of the timelike eigenvector (which 
is e''-^' outside the BH, e^^-' between the outer and inner horizon, and e''-'^' inside the inner horizon), while Pi, Pj, 
and P3, corresponding to the principal pressures in the three spacelike directions, are the three eigenvalues of the 
spacelike eigenvectors. Fig. [I] shows the violation of the weak energy condition for a BH with a/M = 0.6 and 9 = 0.3. 

The second class of transformations (complexification of type-II) includes all the other options. Now m is a function 
of both r and 6*, except in the Kerr limit (7 = 0. The common feature of all these solutions is that there is no global 
transformation Q to write the new metric in the Kerr form in Boyer-Lindquist coordinates. That can happen because 
we are not in the vacuum (assuming the validity of Einstein's equations and introducing some form of exotic matter 
like a non-linear electrodynamics field) and the stress-energy tensor is not the one of a Maxwell electromagnetic 
field, so the hypothesis (m) of Theorem 7.1.1 of Ref. TF* is not satisfied. However, we can still check the absence of 
singularities and the validity of the weak energy conditions, as they are both local properties. The metric is now more 
complicated. We have just considered the case a = —2 and /? = 0. Like for the case a = /3 = 0, it turns out that 
the spacetime is everywhere regular for g 7^ 0, the limit of the curvature invariants assumes different values at r = 
depending on the way one approaches the origin, and the weak energy condition is not respected for a non-vanishing 
spin a. 



IV. BARDEEN BLACK HOLE 



The most famous regular BH solution is the Bardeen metric |4j. In Schwarzschild coordinates, the line element is 
given by Eq. ([T]) with 

/(r) = 1 - — , m = m(r) = M f ' , (26) 

where M is the mass of the BH and g is its magnetic charge. The non-linear electrodynamics field to obtain this 
metric from Einstein's equations was found in Ref. [51. Even in this case, m(r) may be interpreted as the mass inside 



FIG. 1: Plot of p (top left panel), p + Pi (top right panel), p + P2 (bottom left panel), and p + P3 (bottom right panel) for a 
rotating Hayward BH with a/M = 0.6 and g = 0.3. Here x — cos 9. 



the sphere of radius r. The spacetime is everywhere regular (the analytic expression of the curvature invariants is 
reported in IS]) and the weak energy condition is satisfied |S]. 

With the only requirement to recover the Kerr metric for 5 = 0, all the possible complexifications have the form 



2mr , „, / ^2+aj.-a/2 



which is quite similar to the case of the previous section. Once again, we have two classes of solutions. The type-f 
solution has a = /3 — and corresponds to the trivial complexification in which the mass term m{r) is not modified. 



26j) replacing M, 
o not recover the 



Such a solution can be written in the Kerr form in Boyer-Lindquist coordinates with m{r) in Eq 
and it is of Petrov type D. All the other complexifications lead to solutions of type-II, in which we c 
Boyer-Lindquist form of the metric. Both type-I and type-II BHs seem to be everywhere singularity free (that is true 
for the type-I solution, we have checked it is true even for some type-II, and we guess it is true for all these metrics). 
The curvature invariants always assume different values for r = 0, depending on the way we approach the origin. For 
instance, the type-I solution has 

/ \ 2AM [ \ IUAP [ \ 96M2 
lim lim i? U= - — , lim hm R.^R^'' = ^ , lim hm R.^naR^'"'"' = ^ , (28 

while the limit is zero otherwise. The weak energy condition of all these solutions is violated when a 7^ 0. 



V. REVISING THE NEWMAN-JANIS ALGORITHM 



On the base of the results presented in the previous sections, it seems like the Newman- Janis algorithm preserves 
the singularity free property of the non-rotating solution, but not the weak energy condition. Moreover, a peculiar 
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feature concerning the limit r — > of the curvature invariants shows up. These two properties are not very appeaUng 
and we may think about the way to avoid them. 



As already pointed out in the previous sections, the quantity m(r) in Eqs. (15) and (26) looks like the mass inside 
the sphere of radius r and it reduces to M, the BH mass, at large radii. In the third step of the Newman- Janis 
algorithm, we introduce the quantity a by hand, see Eq. (|6]), and we then identify such a parameter with the specific 
spin angular momentum of the BH, i.e. a — J/M. In the standard case of Schwarzschild and Kerr BHs, in which the 
mass is concentrated at the origin r = 0, there are no ambiguities. However, for all the regular BHs in the literature 
the mass of the object seems to be smeared over a larger volume, and it is thus questionable that the specific spin is 
independent of r. 

We have tried to apply this idea to the type-I rotating Hayward BH solution, as it is the one with the simplest 
analytical form. This passage is clearly arbitrary, as we do not see any natural choice for a(r) and we do not know 
when the spin parameter should be promoted to the rank of function. We have thus considered the simplest option, 
in which the rotating BH solution looks like the Kerr metric 

ds^ ^ n — j dt^ H dtd(l) - —dr^ - lldO^ ~ sin^ 9 ir^ + a'^ ^ j dcfp , (29) 



with m given in Eq. ( 15 ) and a' given by 



9'^ 



(30) 



where a = a'(r — > 00) = J/M is the specific spin angular momentum at large radii and g' is a new constant with a 
role similar to g. This choice of a' is just an example. The new spacetime has the curvature invariants at the origin 
assuming the same values of the non-rotating solution, independently of the way one approaches the point r — Q. 
The weak energy condition is not really satisfied, but the violation can be very small, depending on the value of 5', 
as shown in Figs. [2] and |3] for a rotating Hayward BH with a/M — 0.6, g — 0.3, and g' = 2.0. 

VI. SUMMARY AND CONCLUSIONS 

In this letter, we have applied the Newman-Janis algorithm to the regular Hayward and Bardeen BH metrics. 
In both cases, we have obtained a family of rotating BH spacetimes. For the choice of the complexification in the 
Newman-Janis prescription, we have only required to recover the Kerr solution in the limit of vanishing g. The 
remaining freedom is not an ambiguity of the procedure, but it reflects the fact there is an infinite number of possible 
configurations of the matter field (if we see these metrics as solutions of Einstein's equations with an exotic matter 
field rather than solutions of modified Einstein's equations). Our results can also be extended to other BH solutions. 
Generally speaking, there exists a privileged (and trivial, as it alters the 1/r term only) complexification (type-I 
solution): in this special case, the rotating solution looks like the Kerr metric in Boyer-Lindquist coordinates, with 
the mass M promoted to a function m(r) of the r coordinate only. Independently of the exact form of to, the spacetime 
is of Petrov type D and the motion of a free particle is characterized by the existence of the Carter constant. For all 
the other complexifications (type-II solutions), the new metric can be written in the Kerr form in null coordinates, 
with M replaced by some function rh{r,9), whose exact expression depends on the choice of the complexification: 
here there is no global transformation to recover a solution with a unique non- vanishing off-diagonal coefficient. All 
the rotating solutions generated by the Newman-Janis algorithm from the Hayward and Bardeen metrics seem to 
be singularity free. However, the weak energy condition, satisfied for a = 0, is now violated. These solutions have 
also a questionable behavior of the curvature invariant at the origin. We have thus outlined the possibility that the 
Newman-Janis algorithm should be generalized by introducing a spin parameter function of the radial coordinate. 
Such a possibility sounds quite natural, as all the regular BH solutions have a mass function that may be seen as the 
mass inside the sphere of radius r. Such a modification of the Newman-Janis algorithm avoids the suspicious behavior 
of the curvature invariants at the origin and makes the rotating solution violate the weak energy condition only by a 



small amount. We cannot really exclude that with a proper choice of a' in Eq. (30) the violation of the weak energy 
condition can be completely avoided, but all our attempts to do it failed. 
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